Two-stage melting in systems of strongly interacting Rydberg atoms 
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We analyze the ground state properties of a one-dimensional cold atomic system in a lattice, where 
Rydberg excitations are created by an external laser drive. In the classical limit, the ground state 
is characterized by a complete devil's staircase for the commensurate solid structures of Rydberg 
excitations. Using perturbation theory and a mapping onto an effective low energy Hamiltonian, we 
find a transition of these commensurate solids into a floating solid with algebraic correlations. For 
stronger quantum fluctuations the floating solid eventually melts within a second quantum phase 
transition and the ground state becomes paramagnetic. 
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The quantum melting of solids is a paradigm for the 
manifestation of quantum fluctuations at zero temper- 
ature. This quantum phase transition is driven by the 
competition between the interaction energy giving rise 
to a crystalline structure with localized particles and the 
kinetic energy preferring a delocalization of the parti- 
cles. The qualitative behavior of the transition can be 
described by the Lindemann criterion, which states that 
the solid melts if the fluctuations around the mean posi- 
tion reach a certain fraction of the lattice spacing [l| . The 
zero temperature melting of solids with purely repulsive 
interactions has been studied in detail for Wigner crystals 
[31, and has recently attracted large attention for dipo- 
lar system realized with polar molecules Sj and Rydberg 
gases y-|6| . In this Letter we study the quantum melting 
of crystalline phases in driven Rydberg systems. 

Rapid experimental progress motivated by coherent 
applications like quantum computing [7| and quantum 
simulation [8) has pushed the field of Rydberg atoms 
from single-particle physics into an area where strong 
many-body effects are important J9l4l4l|. Within these 
driven strongly interacting systems, crystalline many- 
body ground states have recently been proposed [4|-l6|, 
where an ordered arrangement of Rydberg excitations 
spontaneously appears in a translationally invariant con- 
figuration of frozen atoms. The absence of a kinetic term 
in the microscopic Hamiltonian describing these systems 
requires a deeper analysis to fully understand the nature 
of the quantum melting of these crystalline structures. 
Here, Rydberg atoms in one-dimensional (ID) lattices 
|15| are particularly suitable to investigate the proper- 
ties of Rydberg crystallization. 

In this Letter, we develop an effective theory allowing 
to study the quantum melting of crystalline phases of 
strongly interacting and driven Rydberg atoms in ID lat- 
tices and determine the complete phase diagram. In con- 
trast to the conventional melting of crystalline solids with 
the quantum fluctuations driven by the kinetic energy, 
here, the quantum fluctuations derive from the compe- 
tition between interaction energy and the external laser 
drive, coupling the atomic ground state to the excited 
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FIG. 1: Ground state phase diagram: We find commensu- 
rate solids with filling / = p/q describing a complete devil's 
staircase for f2 = 0. The lobes with p > 1 are not visible 
within this scale. For increasing Rabi frequency fl, a quan- 
tum phase transition takes place into a gapless floating solid 
(FS) with algebraic correlations for the solid structure and 
in general incommensurate density, and eventually a second 
melting transition occurs into a gapped paramagnet (PM). 
This second transition line satisfles A ~ Q,^^'^'^. 



Rydberg state. Deriving an effective low energy descrip- 
tion, we find a two-stage melting: first from a commen- 
surate solid with true long-range order to a floating solid 
with quasi long-range order, and finally to a paramag- 
netic phase, see Fig. [1] The latter transition from the 
floating solid to the paramagnet is in agreement with 
predictions from a universal scaling function [lg| . 

We start with the microscopic Hamiltonian describing 
the driven Rydberg system: Rydberg excitations are cre- 
ated by two-photon processes via an intermediate atomic 
level that can be adiabatically eliminated. Consequently, 
the internal structure of the atoms can be treated as a 
single hyperflne ground state \g) and a single Rydberg 
level |e), giving rise to a spin 1/2 description. The cou- 
pling between the two states by an external laser field is 
given by the Rabi frequency i} and the detuning A. The 
interactions between the Rydberg states are described 



by a repulsive van der Waals interaction with a Ce coef- 
ficient. Here, we focus on a ID setup, where the atoms 
are trapped in a lattice with lattice spacing a. Using the 
Pauli matrices CTq , the Hamiltonian in the rotating frame 
after applying the rotating wave approximation reads 
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with Pee = (l+cri )/2 [17|. Incoherent processes such as 
spontaneous emission are much slower than typical exper- 
imental timescales [7| and can be ignored. In addition, it 
is important to note that ID lattices suitable for Rydberg 
excitations can be realized using deep optical lattices jl^j , 
optical lOl lll| and magnetic [l9| microtraps, or micro- 
fabricated arrays of thermal vapor cells ^20] . We do not 
require to have one single atom per lattice site; instead 
it is possible to work with collective degrees of freedom 
where a single Rydberg excitation is shared among all 
atoms of one lattice site. In the following we focus on 
the situation with a single atom per lattice site. 

The derivation of the phase diagram from the micro- 
scopic Hamiltonian follows in three steps: (i) we start in 
the classical regime for fi = 0, and derive the devil's stair- 
case structure of the commensurate crystal, (ii) We study 
the influence of the driven dynamics with 17 ^ within 
perturbation theory and derive the effective low energy 
Hamiltonian. (iii) For increasing number of defects, we 
study the effective Hamiltonian within mean-field theory 
and establish the transition into the paramagnetic phase. 

For fl = there are no quantum fluctuations, i.e., the 
Hamiltonian is purely classical. The ground state for 
A > follows a complete devil's staircase of crystalline 
configurations with different commensurate lattice spac- 
ings [2l|, I24I . The stability of each crystalline phase with 
a rational filling factor / = p/q is determined by the 
vanishing of the single particle energy gap, i.e., when it 
is energetically favorable to insert or remove one Rydberg 
excitation. For each / we can calculate the detuning Aq 
in the center of the lobe and its width A^ , which leads 
in the strongly interacting regime with / ^ 1 to 
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Note, that the lobes with with p > \ are extremely small, 
and we will restrict the analysis in the following on the 
stability of the lobes with p = 1, i.e., / = l/q. 

In the following, Xi denotes the position of the i^^ Ryd- 
berg excitation in the lattice. Then, the commensurate 
solid satisfies Xi+i — Xi — q. Within this notation, the 
lowest energy excitations are characterized that the dis- 
tance between two neighboring Rydberg excitations is re- 
duced (increased) by a lattice site, i.e., Xi+i—Xi — q^l; 
the two types of excitations are denoted as particle (hole) 
excitations, respectively. Note, that these elementary ex- 
citations describe fractional spin excitations f21| , as the 
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FIG. 2: Processes in second order perturbation theory, (a) 
Virtual annihilation of a Rydberg excitation, (b) Virtual cre- 
ation of a Rydberg excitation, (c) Hopping of a crystal defect. 



addition (removing) of a Rydberg excitation gives rise to 
q = 1// elementary excitations. Their excitation ener- 
gies in the classical regime are Ep^h — /^ (Au,/2=F(5A) 
for particles and holes, respectively. Here, JA = A — Aq 
is the deviation of the detuning from the center of the 
lobe. In addition, two excitations of equal type have a 
repulsive interaction, while two excitations of different 
type are attractive. 

Next, we focus on the regime with a finite drive 51 7^ 0, 
where the system acquires a highly non-trivial dynam- 
ics. Within this regime, the melting of the commensu- 
rate crystals appears via the nucleation and subsequent 
condensation of defects. For a large detuning A ^ fl, 
the system has a well-defined energy gap and we can de- 
rive an effective low-energy Hamiltonian for the defects 
within perturbation theory in il/A. We first restrict our 
treatment to three classes of states: (i) the crystalline 
ground state |c), (ii) the states \pi) with a particle-like 
defect between the Rydberg excitation Xi and x^+i, i.e., 
Xi^i—Xi = 9 — 1, and (iii) the analogous states \hi) for 
hole-like defects. In second order perturbation theory, 
the effective Hamiltonian contains diagonal terms, shown 
in Fig. [2l^a), and Fig. [2jb), while the particle (hole) de- 
fects also acquire an off-diagonal term corresponding to 
a hopping of the defect, see Fig. [DJc). Adding a constant 
energy such that Ec^O the diagonal terms take the form 
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Here, E^""' and E- ' correspond to the energies of the 
virtual levels depicted in Fig.I^fa-b), which take the form 
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Consequently, the diagonal terms provide an additional 
shift in the excitation energy for the particles and holes. 



which in the limit / ^ 1 reduces to (a — p, h) 



AE^ = {a,\H^e\a,) - {c\H^s\c) -E^=I, 
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with the dimensionless quantities Ip k, — //j k, 0.090. The 
ofF-diagonal defect hopping for holes reduces to 
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where the first process annihilates a Rydberg excitation 
at position Xi with the subsequent creation of the Ryd- 
berg excitation at position Xi ± 1, while the second term 
reverses the order of creation and annihilation. Finally, 
the corresponding term Jp for the hopping of particle 
defects is obtained similarly, and the evaluation of the 
second order processes provides Jp = Kphfl^/A and 
Jh — Khhfl'^/A with Kp w Kh ~ —7/5. We can now 
derive the shape of the lobes for the commensurate solid, 
see Fig. [1] the nucleation of particle like defects destroy- 
ing the commensurate solid takes place at the vanish- 
ing of the excitation gap, i.e., Ep + AEp — 2Jp = 
and determines the upper boundary of the lobes, while 
the nucleation of hole like defects becomes preferable for 
Eh + AEh — 2Jfi — 0, which determines the lower bound- 
ary for the lobes. 

The new phase is characterized by a finite density of 
defects. These defects obey a strong on-site repulsion, 
and consequently, close to the phase transition line the 
density of defects is very low. This allows us to describe 
the qualitative behavior of the quantum phase transition 
within an effective spin-1 theory for the defects 

Hcs ^-JJ2 [s^s-,, + h.c] + U {Stf - ^iSt (6) 
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with the three states \m)i describing the presence of a 
particle (to = 1) or hole (to, = —1) defect at site i. The 
XY interaction with J = Jp k^ J^ includes the hopping 
of the defects and the possibility to create particle-hole 
pairs, while the uniaxial anisotropy U = jhA^jl ac- 
counts for the cost in energy to create defects, and the 
chemical potential [i ~ fib A describes the variation in de- 
tuning. Note, that this model neglects additional weak 
particle-hole symmetry breaking terms, longer range in- 
teractions, higher orders in il/A perturbation theory, 
and the possibility to nucleate several defects at the same 
site. However, these additional terms will only change 
the quantitative behavior of the commensurate lobes. 

The effective one-dimensional spin-1 model gives nat- 
urally rise to a transition from the commensurate solid 
with an excitation gap to a phase with algebraic cor- 
relations {SlSf) — (Sf)^ ~ l/|i — jp^ and linear ex- 
citation spectrum. Varying the detuning bA results in 
a commensurate-incommensurate transition with defects 
behaving as free fermions, i.e., if = 1, while the transi- 
tion at the tip of the lobe is described by the Kosterlitz- 



Thouless universality class with if = 2 23|. How- 
ever, the physical quantity describing the properties of 
the Rydberg atoms is given by the spin-spin correlation 
{Pee Pee ) rather than the defect correlations {SfS^). 
Consequently, we have to provide a mapping allowing to 
calculate the physical quantity from the effective model. 
The defect density Sf can be expressed in terms of the 
position of the Rydberg atoms as 5f — Xi+i ~ Xi — q. 
Consequently, the total defect number Nk — J2]=o ^j ~ 
Xk—XQ—kq between the Rydberg excitation at xq and Xk 
defines the position of the fc*^ Rydberg excitation. As the 
system is translational invariant, we can assume without 
loss of generality xq — 0. Then, the correlations between 
the Rydberg atoms take the form 
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with n = {jii) the defect density in the effective Hamil- 
tonian. Here, Pk{m) denotes the probability distribu- 
tion of Nk to find TO defects. In the regime where 
the system is described by free fermions, the distribu- 
tion function can be determined efficiently numerically 
at short distances using Monte- Carlo simulations with 
correlated random numbers [24j, see Fig. |3l On the 
other hand, the long distance behavior can be derived 
within Luttinger liquid theory [25|, predicting a Gaus- 
sian distribution with a mean value nk and a variance 
K^ = ^(jVfe — nfc)^) — K\og{k/b)/7r^; here b ^ nn denotes 
a short distance cutoff. As shown in Fig. [31 this Gaussian 
distribution even captures the qualitative behavior on 
short distances. Consequently, we find that the Rydberg- 
Rydberg correlations oscillate with period n+q and decay 
on short distances as l/-\/log[j/(n + q)/b], while a tran- 
sition to an algebraic decay occurs at large distances, i.e.. 
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with {Pee ) — !/("- + q)- Therefore, we find that the 
transition from the commensurate solid takes place into 
a novel phase exhibiting quasi long-range order for the 
solid correlation function. As the density is in general 
incommensurate with the underlying lattice structure we 
denote this novel phase as a floating solid. It is important 
to note, that the crossover from the short distance regime 
to the algebraic decay takes place once the different terms 
in Eq. ([7]) start to overlap, i.e., j ~ exp[{n+q)^/2]. While 
this scale is much larger than the typical experimental 
scale achievable, the short distance 1/ ^ylog j/b' decay is 
a unique characteristic for this phase. 

In the regime of large Rabi frequencies with J7 ^ 
Ce/a^ ^ A, the spin Hamiltonian ([T]) reduces to the Ising 
model in a transverse and a longitudinal field. Then, 
the ground state is described by a gapped paramagnetic 
phase with algebraic decay of the correlations as l/r^ due 
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FIG. 3: Decay of the spin correlatio ns from Monte-Carlo sim- 
ulations (crosses) according to c/ y^lnj/b' (solid line) as pre- 
dicted by Luttinger liquid theory {q = 6, n = 0.25). The 
crossover to the algebraic decay takes place at much larger 
distances \i — j\ ^ 10®. The inset shows the oscillations of the 
correlation function. 

to the long range interaction [26]. Consequently, a sec- 
ond phase transition takes place from the gapless floating 
solid to the paramagnet. 

This second transition can be seen as the breakdown 
of an effective theory describing the system in terms of 
defects, and takes place when the number of defects Ui at 
a site i becomes of the order of the spacing between the 
Rydberg atoms q. Note, that this criterion is equivalent 
to classical Lindemann criterion for the melting of solids 
[]| . The effective Hamiltonian at the commensurate den- 
sity /K = describing the qualitative behavior for higher 
defect densities includes multiple defects at a site, and 
also allows for multiple defect hopping, i.e.. 
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with \n)i describing the state with n defects at site i 
(n < for holes and n > for particles). Using mean 
field theory, we obtain for the fluctuations of the defect 
number (nf ) ~ J^/[/^, and comparing this value with the 
mean spacing between the Rydberg atoms q'^ ~ J'^/U^ 
provides the scaling A ~ ^^'^/^^ for the transition line 
from the floating solid towards the paramagnet. Note, 
the influence of the underlying lattice drops out for this 
second transition, and the scaling agrees with the value 
previously derived using a universal scaling function [16j . 
The complete phase diagram exhibiting three different 
phases is sketched Fig. [1] Note that this phase diagram 
is also in good agreement with recent numerical results 
obtained for the spin Hamiltonian ([T]) at high fillings [(| . 
Experimentally, the commensurate solids are the most 
challenging to realize. They can be observed if the crit- 
ical Rabi frequency ilc at the tip of the lobe is much 
larger than the intrinsic decoherence rates, e.g., by ra- 
diative decay. For ^^Rb atoms with a principle quantum 
number n = 80 and a lattice spacing a = 266 nm, we 



find for the g = 32 phase fi^ = 27r x 1.5 MHz. Note that 
isolated radiative decays will lead to an overall reduction 
in Rydberg atom density, but correlations will be largely 
unaffected. 

Finally, the presence of a floating solid with incom- 
mensurate fillings between the commensurate solid lobes 
and the paramagnetic phase is highly remarkable as the 
underlying lattice structure is not required to stabilize 
this phase. Consequently, we expect that this floating 
solid with quasi long-range order also survives a different 
arrangement of the atoms; especially we expect that this 
phase is also present in the 'frozen' Rydberg regime 27J , 
where the atoms are distributed randomly. 
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